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Abstract 

We study the existence and properties of stationary solution of ARCH-type equation 
r t = Qa t , where are standardized i.i.d. r.v.’s and the conditional variance satisfies an 
AR(1) equation er 2 = Q 2 (a + & j r t-j ) + r ja 2 _ 1 with a Lipschitz function Q(x) and 

real parameters a,"f,bj. The paper extends the model and the results in [5j from the 
case 7 = 0 to the case 0 < 7 < 1. We also obtain a new condition for the existence of 
higher moments of r* which does not include the Rosenthal constant. In the particular 
case when Q is the square root of a quadratic polynomial, we prove that r t can exhibit a 
leverage effect and long memory. We also present simulated trajectories and histograms 
of marginal density of at for different values of 7 . 
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1 Introduction 


Doukhan et al. {5j discussed the existence of stationary solution of conditionally het- 
eroscedastic equation 


n = Ct.crt, of = Q 2 (a + ^ b j r t -3 ) > 

3 =1 


( 1 . 1 ) 


where {Ct} are standardized i.i.d. r.v.’s, a,bj are real parameters and Q(x ) is a Lipschitz 


function of real variable x E M. Probably, the most important case of (1.1) is 

Q(x) = \/c 2 + x 2 , 


( 1 . 2 ) 


where c > 0 is a parameter. The model 0 -( |1.2[ ) includes the classical Asymmetric 
ARCH(l) of Engle [7] and the Linear ARCH (LARCH) model of Robinson [ 15] : 


n = Ctcrt, cr t = a + /3^2j d 1 r t - j . 

3= 1 


(1.3) 
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[9] proved that the squared stationary solution {rf} of the LARCH model in (1.3) with bj 
decaying as j d ~ l , 0 < d < 1/2 may have long memory autocorrelations. The leverage effect 
in the LARCH model was discussed in detail in m- Doukhan et al. [5] extended the above 


properties of the LARCH model (long memory and leverage) to the model in (1.1)-{1.2) 
with c > 0 or strictly positive volatility. 

The present paper extends the results of j5j to a more general class of volatility forms: 


n = Ct°t, 


o\ 


t = Q' Z [a + J 2 bjrt-jj + 7°f- i: 

3 = 1 


(1.4) 


where {Ct}> a, bj , Q(x ) are as in (1.1) and 0 < 7 < 1 is a parameter. The inclusion of lagged 
j in (1.4) helps to reduce very sharp peaks and clustering of volatility which occur in 


trajectory of 0-0 near the threshhold c > 0 (see Fig. 1). The generalization from 
0 to ( |1.4[ ) is similar to that from ARCH to GARCH models, see 0, 0, particularly, 


i4| with Q(x ) of (1.2) and bj = 0, j > 2 reduces to the Asymmetric GARCH(1,1) of Engle 


0 . 


Let us describe the main results of this paper. Sec. 2 (Theorems [4] and [5]) obtain sufficient 


conditions for the existence of stationary solution of (1.4) with E|rt| p < 00 and 7 £ [0,1). 
Theorem |4] extends the corresponding result in (0 , Thm. 4) from 7 = 0 to 7 > 0. Theorem[5] 
is new even in the case 7 = 0 by providing an explicit sufficient condition ( 2.24[ ) for higher- 
order even moments (p = 4,6 ,...) which does not involve the absolute constant in the 


Burkholder-Rosenthal inequality (2.11). Condition (2.24) coincides with the corresponding 


moment condition for the LARCH model and is important for statistical applications, see 
Remark [ 2 J The remaining sec. 3-5 deal exclusively with the case of quadratic Q 2 in ( |1.2| ), 
referred to as the Generalized Quadratic ARCH (GQARCH) model in the sequel. Theorem[6] 
(sec. 3) obtains long memory properties of the squared process {rf} of the GQARCH model 


with 7 £ (0,1) and coefficients bj decaying regularly as bj ~ /3j d , j —> 00, 0 < d < 1/2. 
Similar properties were established in [5] for the GQARCH model with 7 = 0 and for 


the LARCH model (1.3) in [9j, [10J. Quasi-maximum likelihood estimation for parametric 
GQARCH model with long memory was recently studied in [13]. See the review paper m 
and recent work [12] on long memory ARCH modeling. Sec. 4 extends to the GQARCH 
model the leverage effect discussed in [S] and m- Sec. 5 presents some simulations and 
volatility profiles for the LARCH and GQARCH models with parameters estimated from 


real data. A general impression from our results is that the GQARCH modification in (1.4) 


(1.2) of the QARCH model in [5] allows for a more realistic volatility modeling as compared 
to the LARCH and QARCH models, at the same time preserving the long memory and the 
leverage properties of the above mentioned models. 


2 Stationary solution 

Denote \/j,\ p := E|£o| p (p > 0), p, p := E^q (p = 1 , 2, ... ) and let 

:= ^ ^ bt-sTs- 


(2.5) 


S<t 
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Since 0 < 7 < 1, equations (1.4) yield 


of = yf 7 l QT‘{a + X t -t) and r t = Ct 


1=0 


\ 


J2YQ 2 (a + X t -t 


( 2 . 6 ) 


e=o 


In other words, stationary solution of (1.4), or 


n = Ct 


\ 


Y YQ 2 (a + Y b i r t-e-j) 

3 =1 


(2.7) 


£=0 


can be defined via (2.5), or stationary solution of 


x t ■= Y b *-^ 


s<t 




Y^Q^ + Xs-i), 


( 2 . 8 ) 


e=o 


and vice versa. 

In Theorem |4] below, we assume that Q in (2.6) is a Lipschitz function, i.e., there exists 
Lipg > 0 such that 

\Q(x) ~ Q(y)\ < Lip g |x- y\, x,y£R. (2.9) 


Note (2.9) implies the bound 


Q 2 (x) <c{ + cfx 2 , 


x G 


( 2 . 10 ) 


where ci > 0, 02 > Lipg and C2 can be chosen arbitrarily close to Lipg. 

Let us give some formal definitions. Let J~t = <r(£ s , s < t),t G Z be the sigma-field gener¬ 
ated by Cs,s < t. A random process {ut,t G Z} is called adapted (respectively, predictable) 
if ut is J^-measurable for each t G Z (respectively, ut is Tt -\-measurable for each t G Z). 

Definition 1 Let p > 0 be arbitrary. 


(i) By L p -solution of (2.6) or/and (2.7) we mean an adapted process {rt,t G Z} with 
E|rt| p < 00 such that for any t G Z the series X t = YC°jL\ b j r t-j converges in L p , the series 
of = YQ 2 ( a + X t _() converges in L p / 2 and (2.7) holds. 


(ii) By L p -solution of (2.8) we mean a predictable process {X t ,t G Z} with E|Ah| p < 00 


such that for any t € Z the series of = YleL o7 £( 3 2 ( a + X t -f) converges in L p / 2 , the series 


YC s <t b t- sCs&s converges in L p and ( |2.8[ ) holds. 
Define 

"Er=ii^r 0 < p < 2, 


Bp ■■= 


T,jLi b2 j) p/2 , p> 2, 


B P:1 : — 


B p /( l- 7 p/2 ), 0 < p < 2, 

B P /( l- 7 ) p/2 , P> 2- 


Note B p = Bpfl. As in |5], we use the following moment inequality, see 0, m, m- 
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Proposition 2 Let {Yj, j > 1} be a sequence of r.v.’s such that E|Y ? -| P < oo for some 
p > 0. If p > 1 we additionally assume that {Yj} is a martingale difference sequence: 
E[lj|Yi,..., ij_i] = 0, j = 2,3,.... Then there exists a constant K p > 1 depending only 
on p and such that 


E IXAI P £ K p 

3 = 1 


ET= 0 < p < 2, 

(E J “i(E|iS-l p ) 2/p ) p/2 , p > 2 - 


( 2 . 11 ) 


Proposition [3] says that equations 


and (2.8) are equivalent in the sense that by 


solving one the these equations one readily obtains a solution to the other one. 


Proposition 3 Let Q be a measurable function satisfying (2.10) with some ci,C2 > 0 and 
{Ci} be an i.i.d. sequence with \p\ p = E|Co| p < oo and satisfying E£o = 0 for p > 1. In 
addition, assume B p < oo and 0 < 7 < 1. 

(i) Let {X t } be a stationary L p -solution of (2.8) and let a t := y/YleLo YQ 2 ( a + ^t-e) ■ 


( 2 . 12 ) 


Then {rt = Ct&t} in ( 2 . 6 ) is a stationary L p -solution of \2.7\ and 

Y\n\ p < C(l + E\X t \ p ). 


Moreover, for p > l, {rt,IFt,t 6 Z} is a martingale difference sequence with 


E[r t \Xt-i\ = 0, E[|r t H F t -i] = 


I p a t 


(2.13) 


(ii) Let {r t } be a stationary L p -solution of (2.7). Then {Xt} in (2.5) is a stationary L p - 
solution of ( 2 . 8 ) such that 


E\X t \ p < CE\r t \ p . 


Moreover, for p > 2 


E [X t Xo] = Erv'^Tbt+sbs, t = 0 , 1 , 


S =1 


Proof, (i) First, let 0 < p < 2. Then E \a t \ p = E\a 2 t \ p > 2 < YT=o \l p,2 \ l E\Q{a + X t _f,\ p < 00. 
Hence, using (2.10), the fact that {Xt} is predictable and \Q(a + X t -e)\ p < |c^ + c|(a + 
X t -z)\ p/2 < C( 1 + \a + X t _ e \ p ) < C( 1 + \X t _(\P) we obtain 

OO 

EH* = \p\ P E\a t \ p < Cj2h p/2 \ e a + E\X t _ e \ p ) 

e=o 

< C{l + E\X t \ p ) < 00, 

proving (H T2j ) for p < 2. Next, let p > 2. Then E|<r t |P < (E^ 0 YE 2 / p \Q(a + X t )\ p ) p / 2 < 
CE|Q(a + Xt)\ p by stationarity and Minkowski’s inequality and hence ( |2.12[ ) follows using 
the same argument as above. Clearly, for p > 1 {rt = Ct&t} is a martingale difference 
sequence and satisfies ( |2.13| ). Then, the convergence in L p of the series in (2.5) follows from 
(2.12) and Proposition [2] 


E| 

3 =1 


-3l P < Cl' 


T,T=i\W> 


X 7 =i® 


2 \p/ 2 


= CB P < 00 . 
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0 <p < 2 
p > 2 




















In particular, C,t^YT= o 7 £ Q 2 (« + Y, s <t b t-e- s r s ) = WXXo 7 £ <2 2 0 + ^-t) = n by the 
definition of rt- Hence, {r*} is a L p -solution of (2.7). Stationarity of {ry} follows from 
stationarity of {X f }. 

(ii) Since {r t } is a L p -solution of ([27]), so r t = ( t a t = Ct\/XXo YQ 2 i a + x t-i) with X t 


defined in (2.5) and {Xt} satisfy (2.5), where the series converges in L p . The rest follows 
as in 0, proof of Prop.3. □ 


Remark 1 Let p> 2 and \p\ p < oo, then by inequality (2.11), {r t } being a stationary L p - 
solution of (2.6) is equivalent to {rt} being a stationary Resolution of (|2.6[) with E|ro| p < oo. 


Similarly, if Q and {(t} satisfy the conditions of Proposition [3] and p > 2, then {Xt} being a 


stationary Resolution of (2.5) is equivalent to {Xt} being a stationary Resolution of (2.5) 
with E|Xo| p < oo. See also ([5j, Remark 1). 

Theorem 4 Let {Ct} satisfy the conditions of Proposition [3] and Q satisfy the Lipschitz 


condition in (2.9). 


(i) Let p > 0 and 


Kp^ p |At|p/ p Li P g B}J P < !, 


(2.14) 


where K p is the absolute constant from the moment inequality in (2.11). Then there exists 


a unique stationary L p -solution {X t } of (2.8) and 

C{p,Q)\p\ p B p 


E\X t \ p < 


1 - K p \p\ p Up p Q B P!l ‘ 


(2.15) 


where C(p,Q ) < oo depends only on p and ci,C2 in (2.10) 


(ii) Assume, in addition, that Q 2 (x) = c 2 + c\x 2 , where Ci > 0, i = 1, 2, and p 2 = E£q = 1. 
Then c 2 l?2,7 < 1 is a necessary and sufficient condition for the existence of a stationary 


L 2 -solution {X t } of (2.8) with a/0. 


Proof, (i) We follow the proof of Theorem 4 in [5|. For n 6 N we recurrently define a 


solution of (2.8) with zero initial condition at t < —n as 


X 


( n ) _ 


0, 


t < — n, 


Es=-n h-sCsVs , t> -n, t E Z, 


(2.16) 


where ai^ := \Jy^=q 7 l Q 2 (a + X^). Let us show that {X [ n ^} converges in L p to a 
stationary Resolution {X f } as n —> oo. 


First, let 0 < p < 2. Let m > n > 0. Then by inequality (2.11) for any t > —m we have 
that 


E|xM_xW| p < KM P { Y, I b t - s \ p E\a^\ p + £ \b t . s \ p E\a^ - a^\ p ] 


—m<s<—n 

=:KM P {S' min + S'Y}- 


—n<s<t 


(2.17) 
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Using | Q(a + x)\ p < C + C3I+, i£i with C3 > C2 > Lip q arbitrarily close to Lip^, see 
proof of Theorem 4, we obtain 


m+s 


S' m ,n < E l&t-3| P E^ /2 ( C + ^ E l X ; 


(™)ip\ 

s-e\ ) 


(2.18) 


—m<s<—n £=0 

Next, using |(£i>o 7+ 2 ) 1/2 “ (E;>o 7+ 2 ) 1/2 | < (£;>(++ - ?/i) 2 ) 1/2 we obtain 

( s+n s+m \ 1/2 

E7'(Q(a + EE ) )-Q(° + E-t)) 2 + E 7V(a + E-]) 

£=0 £=s+n+l ) 


(2.19) 


Hence from the Lipschitz condition in (2.9) we have that 

( s+n 

E7 p£/2 Li P ^E|xE ) -E-tr+ E Y e/2 (c+4E\x^ £ )n 

£=0 £=s+n+l 


s+m 


Combining (2.17) and the above bounds we obtain 

s+m 


E | X (m)_ x (n ) r < £ | 6 ,_ s | p ^ 7 ^/ 2 E|E 3 -E-ll P 


—m<s<t £=0 

s+m 


s+m 


v Z> 2 /2 


+ c y. \+.\ v T.i vin+c E i**—r E v 

—m<s<—n £=0 —n<s<t £=s-\-n +1 

< C'iL p | Ai | p K p +n, 7 + ^lMlpC p E + s> 7 E|E n) -+ (m) l P - ( 2 . 20 ) 

—m<s<t 

where + 7 := Ej=o 7' 7p/2 |& s -ll p , s > 0, n t +n ■= C(l-'f/ 2 )~ 1 K p \n\ p (E l j > t+ n |f>++&fV n , 7 ) -> 
0 (n —> oo). 


Iterating inequality (2.20) as in [5], (6.55) and using K p \p\ p c^ Yl s <t ^t-s, 7 = +’I/+ C 3+V7 < 
1 we obtain lim mjn _ ) . 0O E\xj m> — Xt‘\ p = 0 and hence the existence of X t such that 
lim n ^ oc E|+"' i — X t \ p = 0 and satisfying the bound in (2.15). 


Next, consider the case p > 2. Let m > n > 0. Then by inequality (2.11) for any t > —m 
we have that 


EiE m) -+ (n) r < K p \p\ P 


E ++ 2 /p ki m) i p + E b lsV 2/p \ 


a (m) _ a (n)\p 


—m<s<—n 


—n<s<t 


p/2 


■ Llp|lt|p (i7 mn + R mr ^) 


p/2 


( 2 . 21 ) 


Similarly to (2.18), 


—n— 1 m+s —n— 1 m+s 

l™)l 


R'm,n < E b ts E 7'E 2 +Q(« + X™) \ P < E b t-s E ^ + + c 1e 2/p iE-]i p ) 


s=—m £=0 


s=—m £=0 
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and, using (2.19), 


s+n s+ra 

K,n < E 6 L E2/ lE^ i (Q(a + xS)-Q(a + X^ e )) 2 + E 7V(a + E-i) 


—n<s<t 


£=0 


s+n 


-£=s+n+l 
s+ra 


p/2 


< E bl s ( y Y^^ 2/ nQ(a + X^)-Q(a + X^ i )\P+ E 7 'E 2 ^|Q(a + X^) 

—n<s<t £=0 £=s+n+1 


s+n 


s+ra 


< 


£ i '?-,(Lip^E7E 2/, ’|X<” ) -X<: ) ,|>’+ /(C + C 2 E 2 /»|x‘”>PO) 

n<s<t £=0 £=s+n+l 


—n<s<t 
Consequently, 

E 2 / P | X M _ X H|P < Kt+n + /^/ P |^| 2 / p c 2 £ 6 ^^ E 2 / p | x ( m )) _ x ( n )| P ) 


where **+„ := C(1 - 7) 1 Kp /p \p,\l /p {J2 j>t+n + bi +n7 ) 


—m<s<t 
k2 , l 2 


0 (n —> 00). By iterating 


the last displayed equation and using Kp^ p \jj,\p^ p c^ JX =1 & 2 7 = 7Cp^ p |//|p' /p c|i?2/(l — 7) < 


1 we obtain linim^^oc E 2 / p |X/ m — Xj | p = 0 and hence the existence of Xf such that 
limn^oo E|X/ ' — X t \ p = 0 and satisfying the bound in (2.15). The rest of the proof of part 

(i) is similar as in [5j, proof of Theorem 4, and we omit the details. 

(ii) Note that Q(x ) = y/c 2 + c^x 2 is a Lipschitz function and satisfies (2.9) with Lipg = c2. 
Hence by K 2 = 1 and part (i), a unique L 2 -solution {X t } of (2.8) under the condition 
c 2 ^2,7 = c 2 -Z?2/(1 — 7) < 1 exists. To show the necessity of the last condition, let {Xt} be 
a stationary T 2 -solution of (2.8). Then 


EX 2 = E tf-s E 7^ e Q 2 (° + x s-t) 


S<t 


£=0 

00 


E^E7'E(c? + c 2 (a + X 2 _,) 


S<t 


£=0 


= (H 2 /(l- 7))(c? + c 2 (a 2 + EX 2 )) > c|(H 2 /(l - 7 ))EX 2 

since a / 0. Hence, ^^/(l — 7 ) < 1 unless EX 2 = 0, or {X t = 0} is a trivial process. 
Clearly, (2.8) admits a trivial solution if and only if 0 = Q{a ) = \/c 2 + c^a 2 = 0, or 
ci = C2 = 0. This proves part (ii) and the theorem. □ 

Remark 2 Theorem [4] extends ([5], Thm. 4) from 7 = 0 to 7 > 0. A major shortcoming of 
Theorem [4] and the above mentioned result in [5] is the presence of the universal constant 
K p whose upper bound given in [15] leads to restrictive conditions on R Pi7 in (2.14) for the 
existence of L p -solution, p > 2. For example, for p = 4 the above mentioned bound in m 
gives 

K A iuBl/{\ - 7 ) 2 < ( 27 . 083 )V 4 5 2 2 /(1 - 7) 2 < 1 ( 2 - 22 ) 

requiring B 2 = SyLi b'j to be very small. Since statistical inference based of ‘observable’ 
squares r 2 ,1 < t < n usually requires the existence of Erf and higher moments of 77 
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(see e.g. m ), the question arises to derive less restrictive conditions for the existence of 
these moments which do not involve the Rosenthal constant K p . This is achieved in the 


subsequent Theorem [5j Particularly, for 7 = 0,Lipg = 1 the sufficient condition (2.24) of 
Theorem [5] for the existence of E r^,p > 2 even becomes 


3 =2 


EIOmEn 3 ' < L 


(2.23) 


k= 1 


Condition (2.23) coincides with the corresponding condition in the LARCH case in (HOI, 
Proposition 3). Moreover, (2.23) and (2.24) apply to more general classes of ARCH models 
in (1.1) and ( |1.4| ) to which the specific Volterra series techniques used in [9], [TOJ are not 
applicable. In the particular case p = 4 condition (2.23) becomes 


6-B2 + 4|/j.3 I |&fc| 3 + A4 Ei 6 *i 4<1 . 

k=1 k=1 

which seems to be much better than condition (2.22) based on Theorem [4j 

Theorem 5 Let {Ct} satisfy the conditions of Proposition [3] and Q satisfy the Lipschitz 
condition in (2.9). 

Let p = 2,4,... be even and 

p 


j2( p j )\vj\ u PQj2\ bk \ J < ( l ~^ p/2 - 

j =2 ' J ' k =1 


(2.24) 


Then there exists a unique stationary L p -solution {A - *} of (2.8). 


Proof. For p = 2, condition (2.24) agrees with LipgA?2,7 < 1 or condition ( 2. 14| ) so we shall 
assume p > 4 in the subsequent proof. In the latter case (2.24) implies Lip< 1 and 
the existence of a stationary Absolution { X j} of (2.8). It suffices to show that the above 
Absolution satisfies EX^ < 00. 

Towards this end similarly as in the proof of Thm[4](i) consider the solution { X with 
initial condition at t < —n as defined in ( |2.16[ ). Let a[ n * := 0 ,t < —n. Since E (xj n ^—Xt) 

0 (n —> 00), by Fatou’s lemma it suffices to show that under condition (2.24) 

E(X t H ) p < C, 

where the constant C < 00 does not depend on t, n. 

Since p is even for any t > —n we have that 


zero 

2 


(2.25) 


E(X, 


( n )\P _ 


t -1 


E E[6 t - sl C sA^-'-bt-SrCsA 1 ?] 


si,...,Sp=—n 
p / \ t —1 


E 

3 =2 


P 


s —1 


E E b t-uC,u.a, 


3 j 

s=—n 


(n) 


u u u 


P-J 


(2.26) 
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Hence using Holder’s inequality: 


|E£V“ J 'I < dWl p \t/c\ p ^-MP\ri\ p < c>\— E|er + --Elr^l, 1 < j < P, c > 0 

L pc p p 

we obtain 


nx^r < E(j)iwi4 E iC,i{^E ( <,("y+i^ E ( g 6,-„ad” ) )’’} 

j =2 x ^ 7 s=—n u=—n 

= J2 Pi,t- s E(aW/c 3 y+ Y p 2> tsV{x%y, 

s=—n s=—n 

where X^J := Y^-n h-uC,u (J u \ c ‘A > Lip Q and where 


S — 1 


(2.27) 


: = E“u) WsWvM 

3 =2 y V/ 


3 ’ 


A.t-s e 


p-J/p 


3 =2 


P VJ 


l&t-al|Pjl4- 


The last expectation in (2.27) can be evaluated similarly to (2.26)-(2.27): 

S— 1 


It—1 


e(*St = E E (4 n) y( E 


j =2 x 7 u=—n 
s —1 


(n) 




p-J 


s —1 


< E /5i,t-uE(4 n) /c 3 ) p + E )92,t- u E(xg) p . 


Proceeding recurrently with the above evaluation results in the inequality: 

t -1 

dn) 


E(E n) r < E 


(2.28) 


where 


t—S — 1 


Pt—s •— Pl,t—S ^ ^ ^ ^ 02,t—Ul ' ' ’ @2,t—lLk 

k =1 s<Wfc<*”<wi<t 


Let fa := Ylu=i Pi,t, i = 1 , 2 , /3 := /3 t . By assumption ( [2.24P , 

P / \ OO 

/ 3 i +/? 2 = Er n 4 En j < (!- 7) p/2 

j=2 7^/ fc=l 

whenever <73 — Lipg > 0 is small enough, and therefore 


P 


OO OO 

- n _^Ap /2 E^btl 1 + E /3 2^ ) 


( 1 _ 7 ) p /2 - (! _ 7 ) p /2 

1 Pi 


k =1 


(1 _ 7 ) p /2 1 - / 3 2 


< 1. 


(2.29) 
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Next, let us estimate the expectation on the r.h.s. of (2.28) in terms of the expectations on 


the l.h.s. Using (2.10) and Minkowski’s inequalities we obtain 


s+n 

E 2/p (^ n) ) p < xy E2/p i<2( a+x Si p 
£=0 
s+n 

< '£-l‘ E2/r tf + & a + X ?-e)Y 2 

£=0 

n-\-s 

< c + clYsi^/^x^y, 

£=0 

where C 3 > C 2 > Lipq and C 3 — Lipq > 0 can be chosen arbitrarily small. Particularly, for 
any fixed T e Z 

sup E 2 /P( a ^y < C 3 sup E 2 / p (x | n ^) p + C. 

—n<s<T (1 Tj —n<s<T 


Substituting the last bound into (2.28) we obtain 


sup 

—n<t<T 


E 2 /P( X H)P < P 2/p sup E 2 / p (X^y + C. 
(1 - 7 ) —n<s<T 


(2.30) 


Relations (2.30) and (2.29) imply 


sup E 2 / p (X^ n) ) p < 


C 


—n<t<T 


1 - 


p 2 /p 

( i - 7 ) 


< c 


proving (2.25) and the theorem, too. 


□ 


Example: asymmetric GARCH(1,1). The asymmetric GARCH(1,1) model of Engle 
[7] corresponds to 

of = c 2 + (a + br t -\) 2 + 'ycr 2 _ 1; (2.31) 


or 


a 2 = 6 + ipr t -1 + anr 2 .! + 8 a 2 _ x 


in the parametrization of (|T5], (5)), with parameters in (2.31), (2.32) related by 

d = c 2 + a 2 , 8 = 7 , ip = 2ab, an = b 2 . 


(2.32) 


(2.33) 


Under the conditions that {(^ = i"t/&t} are standardized i.i.d., a stationary asymmetric 
GARCH(1,1) (or GQARCH(1,1) in the terminology of |18j ) process { 77 } with finite variance 
and a / 0 exists if and only if . 62,7 = b 2 /( 1 — 7 ) < 1 , or 

6 2 + 7 < 1, (2.34) 
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see Thm. [4] (ii). Condition (2.34) agrees with condition an+8 < 1 for covariance stationarity 
in [18]. Under the assumptions that the distribution of C,t is symmetric and = EQ < oo, 
provides a sufficient condition for finiteness of Er 4 together with explicit formula 

P40\0(l + an + 8) + ■0 2 ] 


Erf = 


(2.35) 


t (1 — a^/r 4 — 2an8 — <5 2 )(1 — an — (5) 

The sufficient condition of |18] for Er 4 < oo is + 2a\\8 + 8 2 < 1, which translates to 

/r 4 6 4 + 26 2 7 + y 2 < 1 (2.36) 


in terms of the parameters of ( |2.31 ). Condition (2.36) seems weaker than the sufficient 
condition /Z 4& 4 + 66 2 < (1 — 7) 2 of Theorem [ 5 ] for the existence of L 4 -solution of ( |2.31[ ). 

Following the approach in [5], below we find explicitly the covariance function p(t ) : = 
cov(rg,r 2 ), including the expression in (2.35), for stationary solution of the asymmetric 


GARCH( 1 , 1 ) in ( 2 . 31 ). The approach in [ 5 ] is based on derivation and solution of linear 
equations for moment functions m2 ■= Er 2 , 7773(4) := Er 2 ro and 772.4(4) := Er 2 rg. Assume 
that /i3 = E£q = 0, or Er^ = 0. We can write the following moment equations: 

m 2 = (c 2 Ta 2 )/(1 - b 2 -7), 7773 (0) = 0, 

OO 

7773(1) = ^7^E(c 2 + a 2 + 2abr_i + 6 2 r 2 £)ro = 2abm2, 

1 =0 

OO 

7773(f) = 22 7^ E (c 2 + a 2 + 2abr t -e-i + 6 2 r 2 _^_ 1 )r 0 

1 =0 

7-2 

= 2 a 6 ? 7727 i_1 + 6 2 ^7^7773(4 — t— 1 ), 4 > 2 . 


(2.37) 


4=0 


From equations above one can show by induction that 7773(4) = 2067772(7 + 6 2 )* l ,t > 1. 
Similarly, 

7774(0) = /74E ((c 2 + a 2 ) + 2abr 0 + 6 2 rg + 7 < 7 o ) 2 

= /U 4 ^(c 2 + a 2 ) 2 + (2a6) 2 777 2 + 6 4 m 4 (0) + 2(c 2 + a 2 )(6 2 + 7)777.2 + (26 2 7 + Y 

OO 

7774(4) = '22 7^E(c 2 + a 2 + 2a6ri_£_i + 6 2 r 2 _^_ 1 )rg 
4=0 

OOOO OO 

= ^7 £ (c 2 + a 2 )?77 2 + 6 2 ^7 £ ?77 4 (|4 - £ - 1|) + 2a6^7 £ 777 3 (^ - 4 + 1), 4 > 1. 

4=0 7=0 7=7 

Using 2ab Y7Lt 7^ m 3^ - 4 + 1) = 4a 2 6 2 m 2 ^£7 7^(7 + 6 2 )^ _t = 4a 2 6 2 7T7 2 7V(1 - 7(7 + 6 2 )) 
and p(4) = 7774(4) — we obtain the system of equations 

p(0) = 7774(0) — 7772, 


p(t) = b 2 y^7^p(|4 - l - 1|) + 4a 2 6 2 77727V(l - 7(7 + 6 2 )) 
4=0 
7-2 

= b 2 22 7^(4 — ^ — 1) + Cq 4-1 , 4 > 1, 


(2.38) 


4=0 
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where C := b 2 YlltLi 7 1 p(P) + ( 7774 ( 0 ) — m 2 )b 2 + Aa 2 b 2 m 2 'y /(1 — 7(7 + 6 2 )) is some constant 
independent of t and 


7714 ( 0 ) = 


Pa nn.2 


1 — 6 4 /U4 — (2& 2 7 + 7 2 ) 


((c 2 + cT)(l + 6 2 + 7 ) + (2a6)"^. (2.39) 


Note that the expression above coincides with (2.35) given that the relations in (2.33) hold. 

Since the equation in (2.38) is analogous to (2.37), the solution to (2.38) is p(t) = C (7 + 
b 2 ) t ~ 1 ,t > 1. In order to find C, we combine p(t) = C( r y + b 2 ) t ^ 1 and the expression for C to 
obtain the equation C = C 6 2 7 /(l — 7(7 + 6 2 )) + (7774(0) — m%)b 2 + 4a 2 b 2 m2 r Y / (1 — 7 ( 7 + 6 2 )). 
Now C can be expressed as 

= 2 (777.4(0) - m|)(l - 7(7 + b 2 )) + 4a 2 ?7727 

1 - 7(7 + 2 b 2 ) 

together with (2.39) and p(t) = C (7 + & 2 )* ,£ > 1 giving explicitly the covariances of 

process {r 2 }. 


3 Long memory 


The present section studies long memory properties of the generalized quadratic ARCH 
model in (1.4) corresponding to Q(x) = s/ c 2 + x 2 of ( 1 . 2 ), viz., 


n = C t 


\ 


^ 7 ^(c 2 + (o+ ^ bt-i-sTs) 2 )-, te Z, (3.40) 

e=o s<t-i 

where 0 < 7 < l,a/ 0 ,c are real parameters, {Ct} are standardized i.i.d. r.v.s, with zero 
mean and unit variance, and bj,j > 1 are real numbers satisfying 

bj ~ p^- 1 (3 0 < d < 1/2, /3 > 0), (3.41) 

The main result of this section is Theorem which shows that under some additional 


conditions the squared process {t’ 2 } of (3.40) has similar long memory properties as in case 
of the LARCH model (see [9], Thm. 2.2). Theorem [ 6 ] extends the result in (0, Thm. 10) 
to the case 7 > 0. In Theorem @ and below, 0 < 7 < 1, B 2 = b 2 and B(-, ■) is beta 


function. 


Theorem 6 Let {rt} be a stationary L 2 -solution of (3.40)-(3.41). Assume in addition that 
Pa = E[£q] < oo, and E[ 77 ] < 00 . Then 

cov(t-q, t- 2 ) ~ k \t 2d ~ l , t — > 00 

where k\ := 1 — 2d)Erg. Moreover, 


(3.42) 


[nr] 


n d i/ 2 ^( r 2 _ Er 2 ) _^ D[0jl] K2 W d+{1/2) (r), 


n 00 , 


(3.43) 


t =1 


where W d+ (i/ 2 ) is a fractional Brownian motion with Hurst parameter H = d + (1/2) E 
(1/2,1) and n 2 := n 2 / (d(l + 2d)). 
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To prove Theorem [6j we need the following two facts. 

Lemma 7 (j5j, Lemma 12) For aj > 0, j = 1, 2,..., denote 

-Ak ■— a k + ^ ^ ^ ^ • • • Oii p —i p _ 1 oik—i pl k — 1 , 2 ,... . 

0 <p<k 0<ii<---<i p <k 

Assume that YljLi a j < 1 an d 


aj < cj 7 , (3 c > 0 , 7 > 1 ). 

Then there exists C > 0 such that for any k > 1 

A k < Ck “ 7 . 

Lemma 8 Tssnme that 0 < /3 < 1 and 07 ~ cj -7 (3 7 > 0, c > 0). Then 

t -1 

a t ,/3 :=^2(3 j a t -j ~ 

l=o P 


Proof. It suffices to show that the difference D t := — at /(1 — /5) decays faster than a*, 

in other words, that 


t— 1 OO 

A = ^2/3 j (a t - a t -j) -^2/3 j at-j = oft -1 ). 
l=o j=t 

Clearly, J2 t/2<j<t ^ (a t - a t -j) = 0{ / 3 t/2 ) = o(t" 7 ), Y!jLtP 3a t-j = 0{^) = o(t" 7 ). 
Relation J2o<j<t/2 ( a t ~ &t-j) = o(t -7 ) follows by the dominated convergence theorem 
since sup 0 <j < t /2 |ctt — <%-j|t 7 < C and |cti — at_j|f 7 —> 0 for any hxed j > 0. □ 

Proof of Theorem [ 6 | We use the idea of the proof of Thm. 10 in |5]. Denote 


t -1 

^.7 : = 

1=0 

b\ a 

t— 1 

:= ^7 J &; 
1=0 

A := bi-si's, 

S<t 


Wi 

II 


(3.44) 


By the definition of 77 in ( |3.40[ ) we have have the following decomposition (c.f. 0 , ( 6 . 66 )) 


( r t ~ E r 2 t ) - Y, bt-s,^r 2 s - Er 2 ) = 2 al t)1 + U t + V t , y =: ft, (3.45) 

S<t 

where Xt , 7 is the main term and the ‘remainder terms’ Ut and are given by 

OO 

Ut := (C t 2 -EC 2 K 2 , V ta := (3.46) 

£=0 

Vt := 2 Y b t - Sl bt- S2 r Sl r S2 . (3.47) 

S 2 <S!<t 
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Using the identity V t = ( X 2 — EAT 2 ) — J2 s <t^t-2( r t ~ Er 2 ) the convergence in L 2 of the 
series on the r.h.s. of (3.47) follows as in [5j (6.67). Hence, the series for V ) i7 in (3.46) also 
converges in L 2 . 

Let us prove that 


cov(£ 0 , 6 ) ~ 4a 2 cov(X 0i7 , X t)1 ) ~ 4a 2 A 2 t 


2 \2 + 2d-l 


t —> OO. 


(3.48) 


where Af = /3 2 /(1 — 'y) 2 B(d, 1 — 2d). The second relation in ( 3.48| ) follows from 6 <i7 ~ 
(/3/(l — 7 ))t d ~ 1 ,t —> oo, see Lemma [7J and the fact that W i7 = Yh s <t bt-s^fs is a moving 
average in stationary uncorrelated innovations {r s }. Since {Ut} is also an uncorrelated 
sequence, so cov(£o, Ut) = 0 (t > 1), and the first relation in (3.48|) is a consequence of 


E[U 0 X t „] + E[U 0 Vt„] = o(t 2d ~'), 

E[X 0 , 7 V), 7 ] + E[V 0 „(X t „ + V t „)\ = o^- 1 ). 


(3.49) 

(3.50) 


We have E[U 0 X t J = b t ^E[U 0 r 0 ] = O^” 1 ) = o^" 1 ) and E [U 0 V t ^\ = 2 b t ^D t = 0(t d ~ x ) = 
o{t 2d ~ x ), where \D t \ := \E[U 0 r 0 J2 s <o bt- s r s ]\ < EU^(Er^) 1/2 (E(X) s<0 b t - s r s ) 4 ) 1/2 < C fol¬ 
lows from Rosenthal’s inequality in (2.11) since E (^) s<0 bt-sfs) 4 < AhEr^ X ^ s<0 ^t-s) 
< C. This proves (3.49). The proof of (3.50) is analogous to [5] (6.68)-(6.69) and is 
omitted. 


Next, let us prove (3.42). Recall the definition of b-j in (3.44). From the decomposition 


(3.45) we obtain 


It - Er t = i G 


(3.51) 


i =0 


where (pj r/ > 0,j > 0 are the coefficients of the power series 3> 7 (z) := = 

(f - Y^jLi &j i7 ~ J ) _1 , z £ C, \z\ < l given by <^ 0i7 := 1, 


W .7 : = &? l7 + E £ w 

0 <k<j 0<si < — <Sk<j 


b 2 ■ ■ ■ b 2 b 2 7 > 1 

SI,7 Sfc-Sfe-1,7 J — 


From (3.41) and Lemmas [T] and [ 8 ] we infer that 

77,7 = 0 (t 2d ~ 2 ), t 


oo, 


(3.52) 


in particular, <J> 7 (1) = Ylt^oVt ,7 = (1 — 7 )/(l — 7 — B 2 ) < 00 and the r.h.s. of (3.51) is 
well-defined. Relations (3.51) and (3.52) imply that 


cov(r 2 ,To) = ( Pi ( Pj coy (&-ii£-j) ~ 4> 2 (l)cov(&,£ 0 ), i-> 00 , (3.53) 

i,j =0 


see [5], (6.63). Now, (3.42) follows from (3.53) and (3.48). The invariance principle in (3.43) 
follows similarly as in [5], proof of Thm. 10 from (3.51), (3.48) and n ~ d ~ 1 / 2 Ylt^i -^,7 —hD[o,i] 
A 2 H / rf + (i/ 2 )(' 7 '), A 2 = A 2 /d(l + 2 d), the last fact being a consequence of a general result in 
[I]. Theorem [ 6 ] is proved. □ 
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4 Leverage 


For conditionally heteroscedastic model in (3.40) with E Q = E£ 3 = 0, E£ 2 = 1 consider the 
leverage function ht = cov(cr 2 ,ro) = Er/ro, t > 1. Following [ID] , and [5], we say that {ry} 
in (3.40) has leverage of order k > 1 (denoted by {ry} G £{k)) if 

hj <0, 1 < j < k. 


The study in [5j of leverage for model (3.40) with 7 = 0, viz., 


n = Ct C* +{a+ Ybt-sr s y 


t G Z 


was based on linear equation for leverage function: 

h t = 2 am 2 h t + Y b l h t-i + 26 f Y b i+tK, t> 1, 


0<2<i 


i>0 


where m 2 = Er^. A similar equation (4.54) for leverage function can be derived for model 
(3.40) in the general case 0 < 7 < 1. Namely, using Ery = 0, Eryro = r?? 2 l(s = 0),Er 2 ro = 
0 (s < 0),Eror Sl r S2 = l(si = 0)/i_ S2 (52 < «i) as in [S] we have that 


t -1 


h t = Er t 2 r 0 = ^ 7 ^E[(c 2 + ( 0 + Y b t-i- s r s ) 2 )?'o] 

e=o s<t-e 

t- 1 t- 1 

= Y^( 2 am 2 bt -z + b t-£- s E i r sro}) + 2 Y^ b t-£-s 1 b t-e-s2 E [ r s 1 r S 2 r 0 } 

1=0 s<t—i 1=0 S 2 <Si<t—i 

= 2 am 2 b t , 1 + Y + 2 Y hiWhttf, (4.54) 


0 <i<t 


i> 0 


where &t )7 ,6 2 )7 are dehned in ( ]3.44[ ) and uyy )7 := X^=o Y b t-£ b i+t-£- 


Proposition 9 Let {rt} be a stationary L 2 -solution of (3.40) with E|ro| 3 < 00 , < 00 . 

Assume in addition that L> 2 , 7 < 1/5, //3 = E£q = 0. Then for any fixed k such that 
1 < k < 00 : 

(i) if ab\ < 0, abj < 0 ,j = 2 ,k, then {rt} G £{k) 

(ii) if abi > 0, abj > 0, j = 2,..., k, then hj > 0, for j = 1 ,..., k. 


Proof. Let us prove that 


:£*?) 


2 ', 1/2 


< 


t =1 


2|a|m2l?2^ 2 
( 1— 7 )(1 — 3 S 2>7 )‘ 


(4.55) 


Let \b\t i7 := T^lbt-el- By Minkowski’s inequality, 


iY w lh 7 ) V2 ^ < l^|i,7 S 2 1/2 > 


t-i 


1/2 


> 1/2 


(4.56) 


2=1 


^=0 


2=1 
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and therefore | Ylil 1 ^iW^tn\ < ||/i||.B 2 ^ 2 |&|t i7 . Moreover, (XXi ^, 7 ) ^ < B ^ 2 /(I - 7 ), 


,OC 

(E£i ml) 1/2 < B 2 I/2 /( 1-7) and ( E£i(Eo<i<t < m\B%; = ||fc||-B 2 /(l - 


2\l/2 


7 ). The above inequalities together with (4.54) imply 

< 2\a\m 2 Bl /2 /(l - 7 ) + \\h\\B 2 /(l - 7 ) + 2||/ l ||B 2 1/2 B 2 1/2 /(l - 7 ), 


proving (4.55). 


Using (4.54) and (4.55), the statements (i) and (ii) can be proved by induction on k > 1 
similarly to [5j. Since wyi i7 = &i^i+i and = b\, equation (|4.54[) yields 


hi = 2am 2 6i, 7 + 2 = 26i(am 2 + ^/ij6i+i). 


i>0 


i>0 


(4.57) 


> 1/2 


According to (4.55), the last sum in (4.57) does not exceed — 


2|a|m 2 S 2)7 /(l — 3B 2i7 ) < \a\m 2 provided B 2n < 1/5. Hence, (4.57) implies sgn(/ii) = 
sgn(afei), or the statements (i) and (ii) for k = 1 . 

Let us prove the induction step k — 1 —»• k in (i). Assume first that a > 0,bi < 0,b 2 < 
0, • • • ,bk -1 < 0. Then hi < 0 ,h 2 < 0, • • • ,hk~ 1 < 0 by the inductive assumption. By 


(4.54), 


hk — 2(om 2 6fc i7 + ^ ^ hjWj'k,y) T ^ ' ^j, 7 ^ , fc —ii 


i>0 


0 <i<k 


where X)o<i<fc &i, 7 ^fc-* < 0 and | £T >0 him,k,~/\ < ||/i||H 2 / 2 | 6| fc]7 < am 2 |6|fc !7 according to 


(4.55), (4.56). Since bk a < 0 and |fr|fc , 7 = |^fc, 7 | this implies am 2 bk n + J2i>o ^i w t,kr/ < 0, or 
hk < 0. The remaining cases in (i)-(ii) follow analogously. □ 


5 A simulation study 


As noted in the Introduction, the (asymmetric) GQARCH model of (3.40) and the LARCH 


model of (1.3) have similar long memory and leverage properties and both can be used for 


modelling of financial data with the above properties. The main disadvantage of the latter 
model vs. the former one seems to be the fact the volatility at may take negative values 


and is not separated from below by positive constant c > 0 as in the case of (3.40). The 


standard quasi-maximum likelihood (QMLE) approach to estimation of LARCH parameters 
is inconsistent and other estimation methods were developed in 0, 0, m , US- 


Consistent QMLE estimation for 5-parametric long memory GQARCH model (3.40) with 
c > 0 and bj = was discussed in the recent work m The parametric form bj = (3j (i 1 


of the moving-average coefficients in (3.40) is the same as in Beran and Schiitzner [2j for 
the LARCH model. 

It is of interest to compare QMLE estimates and volatility graphs of the GQARCH 
and LARCH models based on real data. The comparisons are extended to the classical 
GARCH (1,1) model 

n = &t(t, cr t = Ju + arf_ 1 + /3a^_ v (5.58) 
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We consider four data generating processes (DGP): 


(L) 

m 

m 

(G) 


LARCH of (1.3), 


(5.59) 


QARCH of (3.40) with 7 = 0, 


QARCH of ( |3.40 ) with 7 > 0, 


GARCH(1,1) of (5.58), 


with standard normal innovations and bj = (3j d . The first three models (L), (Ql), (Q2) 
have long memory and (G) is short memory. The parameters ( a,/3,d ) = (0.0101, —0.1749, 
0.3520) (L), (a,c,/3,d) = (0.0058,-0.0101,0.2099,0.4648) (Ql), {a,c,/3,d, 7 ) = (0.0020, 
-0.0049,0.2394,0.2393,0.7735) (Q2) and (cu,a,/3) = (0.00001, 0.1306,0.8346) (G) are ob¬ 
tained from real data, consisting of daily returns of GSPC (SP500) from 2010 01 01 till 2015 
01 01 with n = 1257 observations in total, by minimizing the corresponding approximate 
log-likelihood functions. The details of the estimation procedure can be found in m- 


Fig. 1 presents simulated trajectories of at of four DGP in (5.59), corresponding to the 


same innovation sequence. Observe that the variability of volatility decreases from top to 
bottom, (Q2) resembling (G) (GARCH(1,1)) trajectory more closely than (L) and (Ql). 
The graph (Ql) exhibits very sharp peaks and clustering and a tendency to concentrate 
near the lower threshold c outside of high volatility regions. This unrealistic ‘threshold 
effect’ is much less pronounced in (Q2) (and also in the other two DGP), due to presence 
of the autoregressive parameter 7 > 0 which also prevents sharp changes and excessive 
variability of volatility series. The graph (G) has different shape and volatility peaks from 
the remaining three graphs which is probably due to the short memory of GARCH(1,1). 
Fig. 2 illustrates the effect of 7 on the marginal distribution of (Q2): with 7 increasing, 
the distribution becomes less skewed and spreads to the right, indicating a less degree of 
volatility clustering. 
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Figure 1: Trajectory of DGP: From top to bottom: (L), (Ql), (Q2), (G). The dashed line in (Ql) 
and (Q2) indicates the threshold c/i/1 — 7 > 0 in (3.401. 



Figure 2: Smoothed histograms of DGP (Q2) for different values of 7 . 
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